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Abstract: Various differential cross-sections are measured in top-quark pair (tt) events 
produced in proton-proton collisions at a centre-of-mass energy of y/s = 7 TeV at the LHC 
with the ATLAS detector. These differential cross-sections are presented in a data set 
corresponding to an integrated luminosity of 4.6 fb“^. The differential cross-sections are 
presented in terms of kinematic variables, such as momentum, rapidity and invariant mass, 
of a top-quark proxy referred to as the pseudo-top-quark as well as the pseudo-top-quark 
pair system. The dependence of the measurement on theoretical models is minimal. 

The measurements are performed on tt events in the lepton-fjets channel, requiring 
exactly one charged lepton and at least four jets with at least two of them tagged as 
originating from a 6-quark. The hadronic and leptonic pseudo-top-quarks are defined via 
the leptonic or hadronic decay mode of the W boson produced by the top-quark decay 
in events with a single charged lepton. Differential cross-section measurements of the 
pseudo-top-quark variables are compared with several Monte Carlo models that implement 
next-to-leading order or leading-order multi-leg matrix-element calculations. 
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1 Introduction 


The large number of top-quark pair (tF) events produced at the Large Hadron Collider 
(LHC) has allowed the ATLAS [1] and CMS [2] experiments to perform precise inclusive and 
differential top-quark related cross-section measurements. Both experiments have recently 
published measurements of the inclusive tt production cross-section in proton-proton (pp) 
collisions at centre-of-mass energies, y/s, of 7 and 8 TeV [3-14] as well as differential 
cross-section measurements as functions of the top-quark transverse momentum (px) and 
rapidity (y), and of the mass (m^f) and y of the tt system [15-17]. These cross-section 
measurements triggered recent work on Quantum Chromodynamics (QCD) calculations of 
heavy quark production [18-23] . 

Precision measurements of tt production provide the opportunity to conduct tests of 
predictions based on perturbative QCD and gain direct information on the gluon parton 
distribution function (PDF) at large momentum fractions (xbj) of about 0.1-0.5 [24]. Dif¬ 
ferential tt measurements are particularly sensitive to the tt production mechanism in QCD 
in a region of parton momentum fractions at large momentum transfers. Such studies can 
lead to improvements in background predictions for Higgs measurements and searches for 
physics beyond the Standard Model (SM). 

In the SM, a top-quark decays to a IT boson and a 6-quark with a branching fraction 
close to unity. Hence there are three tt signatures that correspond to different decay modes 
of the W bosons. The signal for this study is in the single-lepton channel. It corresponds to 
the case where one W boson decays directly, or via an intermediate r decay, into an electron 
or muon and at least one neutrino, and the other into a pair of quarks. The neutrino(s) 
will escape the detector unseen, leading to missing transverse momentum whose magnitude 
is denoted by 

The model dependence tt differential cross section measurements presented at the level 
of top quarks, i.e. corrected for parton shower effects and hadronisation, has been an ongo¬ 
ing concern. This paper presents differential ti cross-section measurements using a defini¬ 
tion where the variables are constructed from an object that is directly related to detector- 
level observables. This top-quark proxy object is referred to as the pseudo-top-quark (t). 
The goal of presenting measurements using a definition where the variables are constructed 
from reconstructed charged lepton, jet and missing transverse momentum objects, is to 
allow precision tests of QCD in final states with top-quarks, using reconstructed objects 
that avoid large model-dependent extrapolation corrections to the parton-level top-quark 
but remain well correlated with corresponding objects reconstructed from the partons. 

In this approach, corrections applied to data depend less on Monte Carlo (MC) models 
that describe the hard scattering process modelled by matrix elements calculated to a given 
order, the emission of additional partons from the hard scattering process, parton shower 
effects, hadronisation and multiple-parton interactions. In particular, at low top-quark p^', 
the modelling of soft parton emissions can significantly modify the kinematics of the top- 
quark pair. Observables based on stable particles can be unambiguously compared to MC 
generator predictions and therefore provide useful benchmarks to test and further develop 
MC models and to adjust free model parameters. 
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The t object can be evaluated for hadronic or leptonic decays of the top-quark from 
the detector information or analogously from the stable hnal-state particles generated by 
MC simulations. The differential cross-sections are measured as functions of the transverse 
momentum pt(^) and rapidity y{t) of the leptonic (ti) and the hadronic (th) pseudo-top- 
quark as well as the transverse momentum PT(iiih); rapidity y{t\th) and invariant mass 
of the reconstructed tt system (fith)- 

This paper is structured as follows. The definition of t together with the detector-level 
or particle-level objects used to construct t are presented in section 2. The correlation in MC 
simulation between reconstructed t observables and corresponding top-quark observables 
at the parton level is also discussed in this section. Section 3 provides a short overview 
of the ATLAS detector. A description of the different MC samples in the study is found 
in section 4. The data and MC event selection is described in section 5 together with 
the reconstruction of final-state objects. Section 6 covers the treatment and evaluation of 
systematic uncertainties. Comparisons between data and MC simulation, for the yields and 
pseudo-top-quark distributions before unfolding, are presented in section 7. A description 
of the unfolding, the results obtained for data and various MC models and concluding 
remarks are found in sections 8, 9 and 10, respectively. 

2 Measurement definition 

The model dependence of parton-level tt cross section measurements has been an ongoing 
concern for QCD studies in particular in events with top quarks. The use of particle- 
based definitions in cross-section measurements is a standard methodology in high-energy 
physics [25, 26] to reduce the model dependence. Standardised tools exist to compare theory 
predictions for such measurements [27-29]. In tt events, such measurements were published 
for the inclusive cross-section [3, 4] and for differential cross-section measurements as a 
function of the transverse momentum and the rapidity of the final-state leptons and jets [17, 
30-32]. In this paper the concept of the particle-based cross-section dehnition is extended 
to the kinematic properties of the top-quark decay products that when correctly combined 
are closely related to the kinematic properties of the top-quark. An operational definition 
t that is defined from measured observables as described below has been introduced to 
reduce the model dependence of the measurements.^ This operational dehnition should 
result in a good correlation between the t object and the top parton for a given Monte 
Carlo generator. However, it is not necessary that the algorithms lead to the best possible 
correlation to the parton-level kinematics or the best signal to background ratio. 

The differential cross-section measurements in this paper are presented in terms of the 
kinematics of the t object, introduced in section 1. The identihcation of reconstructed 
charged lepton, jet and objects from actual or simulated detector signals is discussed 
in section 5.1. The identihcation of particle-level objects for MC events is discussed in 
section 2.1. The kinematic hducial region for both the reconstructed and particle-level 
objects is dehned in section 2.2. The algorithm used to construct leptonic or hadronic 

^ Discussions related to an operational definition have occurred between experimentalists and theorists 
in the Top Physics LHC working group. 
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t objects from either reconstructed detector-level objects or from particle-level objects is 
described in section 2.3. 

The measurements presented in this paper can be directly compared to MC simulations 
using matrix-element calculations for the hard scattering, interfaced with parton shower and 
hadronisation models. These are referred to as particle-level predictions. For comparisons 
to fixed-order QCD calculations, corrections for the transition from partons to hadrons 
need to be applied. 

2.1 Particle objects 

In the case of MC simulation, objects can be identihed at the particle level. Leptons and jets 
are defined using particles with a mean lifetime r > 3 x 10“^^ s that are directly produced 
in pp interactions or from subsequent decays of particles with a shorter lifetime. The 
lepton dehnition only includes prompt electrons, muons and neutrinos not originating from 
hadron decays as well as electrons, muons and neutrinos from tau decays. The electron 
and muon four-momenta are calculated after the addition of any photon four-momenta, 
not originating from hadron decay. The electron and muon four-momenta are defined 
to include any photons not originating from hadron decays that are found within AR = 
\J{A(j)Y + (Ar/)2 <0.1 with respect to the lepton direction. The direction of the photons is 
required to be within AR = (Aip)^ + (At/)^ = 0.1 with respect to the lepton direction 
The missing transverse momentum vector and its associated azimuthal angle are evaluated 
from the sum of the neutrino four-momenta, where all neutrinos from W boson and r 
decays are included. Jets are dehned by the anti-fcf algorithm [33] with a radius parameter 
of 0.4. The jets include all stable particles except for the selected electrons, muons and 
neutrinos, and the photons associated with these electrons or muons. The presence of one 
or more 6-hadrons with p^ > 5 GeV associated to a jet defines it as a 6-jet. To perform the 
matching between 6-hadrons and jets, the 6-hadron energy is scaled to a negligible value 
and included in the jet clustering (ghost-matching) [34]. 

2.2 Kinematic range of objects 

The cross-section measurement is defined in a kinematic region where the reconstructed 
physics objects have a high reconstruction efficiency (hducial region). The kinematic region 
is chosen such that the kinematic selections of the physics objects reconstructed in the 
detector and of the particle objects are as close as possible. The hducial region is dehned 
in the same way for reconstructed physics and particle objects. However, on detector-level 
some additional selections can be applied. 

Electrons, muons and jets are required to satisfy px > 25 GeV and \r]\ < 2.5. The 
hducial volume is dehned by requiring exactly one muon or electron, four or more jets of 
which at least two are 6-jets, E™'®® > 30 GeV and a W boson transverse mass m'YiW) > 

^ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre 
of the LHC ring, and the j/-axis points upward. Cylindrical coordinates (r, <()) are used in the transverse 
plane, 4> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as g = — lntan(^/2). 
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35 GeV.^ Events are discarded if the electron or muon is within AR = 0.4 of a jet, or two 
jets are within AR = 0.5 of each other. 

2.3 Hadronic and leptonic pseudo-top-quark definition 

The definition of t as a hadronic or leptonic object is determined by the decay of the W 
boson. In the cross-section definition used in this paper, the two highest pT ^-jets are 
assumed to be the 6-jets from the top-quark decay. 


In the case of ti, the leptonically decaying W boson is constructed from the electron 
or muon and the The 6-jet with the smallest angular separation {AR) from the 

electron or muon is then assigned as a decay product of ti. Using the measured W 
boson mass, mw = 80.399 GeV [35], and the components of the missing transverse 
momentum vector (denoted as Px,u,Py,u) associated with the W boson decay neutrino, 
the Pz^u of the neutrino can be constrained: 


T Ejz) {Px,£ T Px,u) iPy,£ T Py,i^) iPz,l T Pz,u) — 


m 


W’ 


where the subscript £ refers to the electron or muon. Neglecting the neutrino mass, 
the Pz^iy of the neutrino is taken from the solution of the resulting quadratic equation: 


where 


and 


Pz,i. = 


—6 ± y/lP' — 4ac 


2a 


a = Ej- pI , 


6 = -2kpz,e c = Ejp\ ,^ - k 


k = 




T {Px,lPx,u T Py,£Py,i') ■ 


If both solutions are real, the solution with the smaller \pz,u\ is chosen. In cases where 
(6^ — 4ac) is less than zero, pz^u is taken as: 


2a- 

Given the value of Pz,u, b is formed from the combination of the charged lepton, 
neutrino and assigned 6-jet. 

In the case of the hadronically decaying W boson is constructed from the remaining 
two highest-pT jets. The is then defined from the hadronically decaying W boson 
candidate and the remaining 6-jet. 


®The W boson transverse mass mT{W) is defined as — cos{(j>^ — 4‘'^)), where I and v refer 

to the charged lepton (e or /r) and the missing transverse momentum vector, respectively. The symbol (j> 
denotes the azimuthal angle of the lepton or missing transverse momentum vector. 
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Once the leptonic and hadronic t are defined, their four-momenta can be evaluated, 
and used in the measurement to define, for example the transverse momentum (pT(ii) 
or PT(ih)) the rapidity (y(ti) or y(th)) and the transverse momentum (pT(tith))) rapidity 
(y(tith)) and invariant mass {m{tith)) of the reconstructed tt system. 

In summary, three different definitions are used in the following sections: 

• A parton-level top-quark is the MC generator-level top-quark selected before it decays 
but after any radiative emissions;'^ 

• A particle-level pseudo-top-quark (hadronic and leptonic) is defined by stable generator- 
level particles within the described acceptance; 

• A detector-level pseudo-top-quark (hadronic and leptonic) is evaluated with the use of 
physics objects reconstructed from detector measurements as discussed in section 5.1. 

Detector-level distributions of t variables are corrected for detector efficiency and reso¬ 
lution effects to allow comparison with distributions of equivalent particle-level t variables 
(see section 8). Therefore any MC model that simulates the final-state particles from pp 
collisions can be compared to these data. This makes existing or future MC model compar¬ 
isons possible, independent of the presence of top-quark partons in the MC event record. 
For comparisons of data to simulation it is not necessary to use the tt kinematic recon¬ 
struction method with the best performance, but it is important that the method is well 
defined and the same definition is applied to data and simulation. 

2.4 Kinematic comparison of the parton-level top-qnark with the psendo-top- 
qnark 

The pseudo-top-quark definition is chosen such that it is closely related to the top-quark 
parton provided by QCD calculations. 

Figure 1 compares the parton-level pT distribution of the top-quark with the particle- 
level pt distribution of the pseudo-top-quark. The two distributions are shown for the same 
integrated luminosity, using the powheg-|-pythia Monte Carlo generator (see section 4). 
The number of events with a reconstructed t is much smaller than the total number of 
generated tt events. This is primarily due to the requirement of a single lepton, four jets 
and the pseudo-top construction efficiency. 

The ratio of the normalised distributions illustrates the shape difference. The larger 
phase space of the parton-level distribution results in a softer px distribution when com¬ 
pared with the pseudo-top-quark distribution where fiducial cuts are applied. In addition, 
with increasing pxj the decay products tend to be more collimated due to their large boost 
resulting in a reduced pseudo-top reconstruction efficiency. Using the same Monte Carlo 
sample, figure 2(a) shows the correlation between the generated parton-level px distribu¬ 
tion and the particle-level hadronic pseudo-top-quark px- About 50% of the events show 

The four-momenta of the top-quarks produced in the hard scattering 2^2 process are modified after 
gluon emission such that an event record with consistent information and conserved energy and momentum 
balance is obtained. Technically, the parton level is defined as status code 155 for herwig and status code 
2 for PYTHIA. 
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Figure 1. Simulated event distribution of the parton-level top-quark px distribution for all events 
(green), and the particle-level pseudo-top-quark px (^h) for events within the fiducial region (blue). 
In both cases the top-quark that decays hadronically is chosen. The distributions are evaluated for 
the same event sample based on powheg-|-pythia at y/s = 7 TeV. The upper hgure is made for 
an arbitrary integrated luminosity. The lower figure shows the ratio of the particle-level th over the 
parton-level top-quark normalised distributions to emphasise the difference in shape between the 
two. 

a strong correlation between the top quark parton and the particle-level th- These are the 
cases where the correct combination of the final state object was found to reconstruct the 
particle-level ih- Despite the good correlation, there remains a significant bin migration in 
the region populated by the majority of top-quarks. These migrations can be affected by a 
change in the modelling of tt production. This makes any measurement extrapolated to the 
parton level model-dependent. However, there is a much stronger correlation between the 
the hadronic particle-level and detector-level pseudo-top-quark px distributions as shown 
in figure 2(b). Hadronic pseudo-top-quark measurements that are presented at the particle 
level are therefore less affected by model-dependent corrections. 

3 ATLAS detector 

The ATLAS detector [1] is a general-purpose detector that covers nearly the entire solid 
angle around one of the pp interaction points of the LHC [36]. It is composed of an inner 
tracking detector (ID), covering a range of \r]\ < 2.5, surrounded by a superconducting 
solenoid that provides a 2 T magnetic field, high-granularity electromagnetic (EM) and 
hadronic sampling calorimeters and a muon spectrometer (MS) that incorporates a system 
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Figure 2. (a) Monte Carlo study using the nominal powheg-I-pythia MC sample showing 

the correlation between the parton-level top-quark px and the particle-level hadronic pseudo-top- 
quark pt and (b) the correlation between the particlc-level hadronic pseudo-top-quark px and the 
hadronic pseudo-top-quark px evaluated from reconstructed objects. In each case the correlation is 
normalised to all the events within a bin on the horizontal axis. 


of three air-core superconducting toroid magnets, each with eight coils. The ID comprises 
a silicon pixel detector, a silicon microstrip detector (SCT), and a transition radiation 
tracker (TRT). The EM barrel calorimeter is composed of a liquid-argon (LAr) active 
medium and lead absorbers. The hadronic calorimeter is constructed from steel absorber 
and scintillator tiles in the central pseudorapidity range of \rj\ < 1.7, whereas the end-cap 
and forward regions are instrumented with LAr calorimeters for both the electromagnetic 
and hadronic energy measurements up to |p| = 4.9. The MS toroid magnets are arranged 
with an eight-fold azimuthal coil symmetry around the calorimeters. Three layers of muon 
spectrometer chambers surround the toroids. High-precision drift tubes and, at small 
radius in the end-cap, region cathode strip chambers provide an independent momentum 
measurement. Resistive plate chambers in the central region and fast thin gap chambers 
in the end-cap region provide a muon trigger. 

Data are selected from inclusive pp interactions using a three-level trigger system. A 
hardware-based first-level trigger is used to initially reduce the trigger rate to approximately 
75 kHz. The detector readout is then available for two stages of software-based (higher- 
level) triggers. In the second level, partial object reconstruction is carried out to improve 
the selection and at the last level, the event filter, a full online event reconstruction is made 
to finalise the event selection. During the 2011 run period, the selected event rate for all 
triggers following the event filter was approximately 300 Hz. 














4 Monte Carlo simulation 


Monte Carlo simulation samples were generated to correct the measurements for detector 
effects. The production of tt events is modelled using the powheg [37-40], MC@NLO [41] 
and ALPGEN [42] generators, powheg and mg@nlo use NLO matrix-element calculations 
interfaced with parton showers and alpgen implements a leading-order (LO) multi-leg 
calculation for up to five additional partons, which are subsequently matched to parton 
showers [43]. A common feature of the samples, unless stated otherwise, is the generation of 
the underlying event and parton shower, which is either performed by pythia (v6.425) [44] 
or by HERWIG (v6.520) [45] together with jimmy (v4.31) [46]. The program tauola [47] 
is used for the decays of r-leptons and photos [48] for photon radiation. 

A POWHEG sample was generated using the version powheg- hvq p4 [39], with the 
CtIOnlo pdf set [49] and the default factorisation and renormalisation scales set to 
where mt is the top-quark mass and pT the top-quark transverse momentum 
as evaluated for the underlying Born configuration (i.e. before radiation). The powheg 
matrix-element calculation was interfaced with pythia, using the “C” variant of the Perugia 
2011 tunes [50] and the corresponding gteqGlI PDF set [51]. This sample is referred to 
as “powheg-|-pythia” and is the benchmark signal sample in this study. 

An additional powheg sample was generated in order to compare the parton shower 
and fragmentation models. It was interfaced with herwig together with jimmy for the un¬ 
derlying event model using the AUET2 tune [52] (“powheg-|-herwig”). To evaluate the im¬ 
portance of the gluon PDF on the final corrected distributions, an alternative powheg sam¬ 
ple was produced with the herapdf15nlo PDF set [53] (“powheg(herapdf)-I-pythia”). 
This PDF set is based on HERA I data together with the inclusion of the precise high-Q^ 
preliminary HERA H data. Simulations using this PDF set are in good agreement with 
early W and Z boson production measurements at the LHC [54]. 

To compare the powheg predictions, alternative samples were generated using MG@NLO 
(v4.01) with the gtIOnlo PDF set and its default renormalisation and factorisation scales: 

= {{PTt +PTf)/^) + mc@nlo was interfaced to herwig and jimmy, using the 
AUET2 tune (“mg@nlo-|-herwig”). 

To assess the modelling of LO matrix-element calculations for additional partons, the 
alpgen [42] generator (v2.13) was used together with the GTEQGlI PDF set and associated 
strong coupling constant, as (mz) = 0.129. The produced processes correspond to the LO 
matrix elements for tt with five inclusive associated partons and tt + bb and tt + cc states, 
at its default renormalisation scale and with the factorisation scale set to ^ m? +p‘^. 

Here the sum runs over heavy quarks and light jets® with mass m and transverse momentum 
Pt- The ALPGEN samples were interfaced with herwig and jimmy, using the MLM parton- 
jet matching scheme [42] with a matching scale of 20 GeV. The exclusive heavy-flavour 
samples were combined with the inclusive samples, after the removal of overlapping events. 
This sample is referred to as “alpgen-|-herwig”. 

^Defined to be a jet comprising gluons and light quarks and subsequently used in the MLM matching 
scheme. 
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Both the NLO matrix-element-based MC models and the LO multi-leg MC models 
have higher-order correction uncertainties that can be estimated in terms of initial-state 
radiation (ISR) and final-state radiation (FSR) variations, alpgen (v2.14) is used to 
generate ti samples with the CTEQ5l PDF set [55], the pythia parton shower and the 
Perugia 2011 tune. Nominal and shifted ISR/FSR samples were produced with an ag 
corresponding to Aqcd = 0.26 GeV as used in the Perugia tune and by modifying the 
renormalisation scale at each local vertex in the matrix element by a factor of 2.0 (0.5) 
relative to the original scale to obtain more (less) radiation. The renormalisation scale is 
varied by the same factor in the matrix-element calculation and pythia, where the radHi 
and radLo Perugia 2011 pythia tunes are used [50]. These samples are referred to as 
ALPGEN-|-pythia(q!s Up) and ALPGEN-|-PYTHlA(as Down), respectively. The selected ag 
values are found to produce variations that are similar to the uncertainty band of cross- 
section measurements for tt events with additional jet activity, as described in ref. [30]. 
The effect of colour reconnection is estimated by generating a powheg-|-pythia sample 
in which no colour reconnection is allowed within PYTHIA, using the noCR Perugia 2011 
tune [50]. 

The total inclusive tt cross-section for pp collisions at y/s = 7 TeV is calculated to 
be ati = 1771J? pb for a top-quark mass mt = 172.5 GeV. This calculation is carried 
out at next-to-next-to-leading order (NNLO) in QCD including resummation of next-to- 
next-to-leading logarithmic (NNLL) soft gluon terms [56-61] with top-| — h (v2.0) [62]. All 
tt MG samples were generated with nit = 172.5 GeV and were normalised to the theoretical 
NNLO-I-NNLL cross-section prediction. 

For the simulation of background processes, samples of W and Z bosons produced in 
association with jets were generated using alpgen (v2.13) with LO matrix elements for up 
to five inclusive associated partons. The GTEQ6l1 PDF set and the herwig parton shower 
were used. In addition to the inclusive jet-flavour processes, separate samples of IF65-|-jets, 
VFcc-|-jets, kPc-l-jets and Zbb+iets matrix-element processes with three additional partons 
were generated and the overlap between them removed. The normalisation of the IT-l-jets 
samples were determined from data as described in section 5.3. The Z-|-jets samples are 
normalised to the cross-section obtained from an NLO QGD calculation with MGFM [63] 
using the mstw2008nlo PDF set [64]. 

A sample of t-channel single top-quark decays was generated using the agermg gener¬ 
ator [65] (v3.8), while MG@NLO was used to generate VFt-channel and s-channel processes. 
Each of these samples is normalised according to an NLO-I-NNLL calculation for the cor¬ 
responding t-channel [66], s-channel [67] and VFt-channel [68] processes. Diboson events 
iWW, WZ, ZZ) were produced using herwig and normalised to the cross-section obtained 
from an NLO QGD calculation with MGFM using the mstw2008nlo PDF set. 

To properly simulate the LHC environment, additional inelastic pp interactions were 
generated with pythia using the AMBTl tune [69] and then overlaid on the hard process. 
The MC events are weighted such that the distribution of the generated mean number of pp 
collisions per bunch crossing ((/r)) matched that of the data-taking period. The particles 
from additional interactions are added before the signal digitisation and reconstruction 
steps of the detector simulation, but are not used within the particle-level measurement 
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defined in section 2. 

The response of the detector to the generated events is determined by a full geant4 [70] 
simulation of the ATLAS detector [71]. This is performed for all samples except for the 
ISR/FSR variations, colour reconnection and powheg+herwig MC samples. For those 
samples a faster simulation which parameterises the ATLAS calorimeter response is used 
instead [71]. 

5 Data sample and event selection 

The data are selected from the full 2011 data-taking period. Events are required to meet 
baseline data quality criteria during stable LHC running periods. These criteria reject 
data with significant detector noise or read-out problems and depend on the trigger con¬ 
ditions and the reconstruction of physics objects. The resulting data set corresponds to 
an integrated luminosity of 4.59 ± 0.08 fb“^ [72]. During this period, the LHC delivered 
instantaneous luminosities that were sufficient to produce several pp interactions within the 
same bunch crossing (in-time pile-up). Interactions in adjacent bunch crossings also influ¬ 
enced the detector and readout signals in the selected bunch crossing. The mean number 
of in-time pile-up interactions, (/i), was measured by averaging over all pp bunch crossings 
in a given luminosity block. The average value of (p) was approximately 5 at the beginning 
of the data-taking period and as high as 18 by the end of the 2011 run. 

5.1 Object reconstruction 

Primary vertices are formed from tracks reconstructed in the ID. The selected primary 
vertex is required to include at least four reconstructed tracks satisfying px > 0.4 GeV and 
to be consistent with the pp beam collision region in the transverse plane. In the cases 
where more than one primary vertex with at least four tracks is reconstructed, the vertex 
with the highest ^ of the associated tracks is chosen and assumed to be associated with 
the hard process. 

Electron candidates are identified as electromagnetic energy deposits (clusters) matched 
to a reconstructed track in the ID [73] . Selected electrons are required to satisfy stringent 
identification criteria. The reconstructed tracks are required to have a minimal number 
of pixel and SCT hits among those expected along the electron trajectory and at least a 
minimum number of high-threshold TRT hits. The longitudinal and lateral shower profiles 
in the calorimeter are required to match those expected for an electron, with a satisfactory 
match between the cluster energy and the reconstructed track momentum. 

To reduce the rate of events with non-prompt and fake lepton signatures from multi¬ 
jet background processes, electrons are required to be isolated within both the calorimeter 
and ID. The calorimeter isolation is defined using a cone of size AR = 0.2 around the 
electron direction. The transverse energy sum of the clusters found in the cone is calculated 
and required to be less than 10% of the electron transverse energy, after excluding the 
calorimeter cells associated with the electron cluster and correcting for leakage from the 
electron cluster. The ID-based isolation is calculated using the summed track-px within 
a AR = 0.3 cone around the electron direction and is required to be less than 10% of 
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the electron track px- Electrons are selected by requiring px > 25 GeV in the range 
\r]\ < 2.47, excluding the barrel/end-cap transition region of 1.37 < \ri\ < 1.52. Electrons 
with px > 15 GeV are used for the object overlap removal discussed later in this section 
and to remove events with two or more leptons as discussed in section 5.2. 

Muon candidates are required to be composed of a reconstructed track in the MS 
combined with a track in the ID [74]. Track quality criteria are used to reduce the non¬ 
prompt and fake muon background from multi-jet processes and to select a sample with 
improved px resolution. Reconstructed tracks are required to have a hit in the innermost 
pixel layer if expected from the track trajectory and at least a minimum number of pixel 
and SCT hits, set below the number of expected hits on a muon trajectory. Muons crossing 
the TRT are required to have a hit pattern consistent with a well-reconstructed track. 

To further reduce the fake muon background, muons are required to be isolated within 
the calorimeter and ID. The calorimeter isolation is determined using using calorimeter 
energy deposits in a AR = 0.2 cone around the direction of the muon and is required to 
be less than 4 GeV. The ID isolation is determined from the summed px of tracks in a 
AR = 0.3 cone around the direction of the muon and is required to be less than 2.5 GeV, 
excluding the px of the muon. The muon channel event selection requires the reconstruction 
of one muon with px > 25 GeV associated with the selected primary vertex. Muons with 
Px > 15 GeV are used to define an additional lepton veto discussed in section 5.2. Both 
types of muons are selected within jpj < 2.5. 

Topological clusters [75] are formed from calorimeter energy deposits. Jets are recon¬ 
structed from these clusters with the anti-fcf algorithm [33] with a radius parameter of 0.4. 
The jets are calibrated using the EMj-JES scheme, where the jet energy scale (JES) is 
derived as a correction of the initial calorimeter calibration set for electromagnetic showers 
[76]. The jet energy is corrected for the effect of additional pp collisions in data and MG 
events. 

To correct for energy losses due to non-compensation in the calorimeter, uninstru¬ 
mented material and detector subsystems in front of the calorimeter, jet energy corrections 
factors are applied that depend on the jet energy and the jet r/ to achieve a calibration 
that matches the energy of stable particle jets in simulated events (excluding neutrinos and 
muons). 

Differences between data and Monte Carlo simulation are evaluated using in situ tech¬ 
niques and are corrected in an additional step [77]. The in situ calibration exploits the px 
balance in events with a Z boson {Z+]ei) or a photon (y-l-jet) and a recoiling jet and dijet 
events. Z-|-jet and y-l-jet data are used to set the JES in the central detector region, while 
Px balancing in dijet events is used to achieve an rj intercalibration of jets in the forward 
region with respect to central jets 

The calibrated jets are required to satisfy px > 25 GeV and be within the range 
IpI < 2.5. Jets associated with large energy deposits from additional pp interactions are 
removed by requiring that the px sum of the reconstructed tracks matched with the jet 
and the selected primary vertex is at least 75% of the total px sum of all tracks associated 
to the jet. This quantity is referred to as the jet vertex fraction (JVE). Jets satisfying 
Px > 50 GeV are always accepted and jets having no associated tracks are also accepted. 
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The MVl algorithm [78] is used to select jets associated with 6-hadron decays. The 
algorithm combines several tagging algorithms into a single neural-network-based discrim¬ 
inant. Jets are identified as 6-jets by using an MVl discriminant tuned to achieve a 70% 
tagging efficiency for jets with px > 20 GeV in simulated tt events. The corresponding re¬ 
jection factor for jets originating from gluons or light quarks is found to be approximately 
130. 

The and its associated azimuthal angle are reconstructed from the vector sum of 

the transverse momenta of the reconstructed objects (electrons, muons, jets) as well as the 
transverse energy deposited in calorimeter cells not associated with these objects, within 
the range |r/| < 4.9. The object classification scheme for the electrons, muons and jets used 
to calculate E™®® is chosen to be the same as the object dehnitions used in this analysis. 
Calorimeter cells not associated with an object are calibrated at the electromagnetic scale 
before being included in the E™^®® calculation. This calibration scheme is similar to the 
one described in ref. [79]. 

Electron and jet objects are reconstructed using separate algorithms that are run in¬ 
dependently. Jets are reconstructed from topological clusters as described above, with no 
distinction made between identified electron and jet energy deposits within the electro¬ 
magnetic and hadronic calorimeters. Jets associated with an electron energy deposit are 
discarded using angular matching. For each electron, the jet with an axis closest to that 
of the electron direction, within AR < 0.2, is discarded. To remove leptons from heavy- 
flavour decays, the lepton is discarded if the lepton is found to be within AR < 0.4 of a 
selected jet axis. 

5.2 Event selection 

Data were collected by requiring either a high-px electron trigger, based on calorimeter 
energy deposits, shower shape and track quality constraints, or a high-px muon trigger that 
included a reconstructed track in the MS matched with a track in the ID. The electron 
trigger px threshold was either 20 GeV or 22 GeV, depending on the data-taking period, 
whereas the muon trigger px threshold remained at 18 GeV for the duration of the entire 
2011 data taking period. 

The selected events are required to contain at least one reconstructed primary vertex. 
A small number of events are rejected that included one or more jets of px > 20 GeV with 
energy that is identihed as being from noise in the calorimeter electronics, from non-pp 
collision background sources or cosmic-ray showers. Events where an identified electron 
and muon share the same reconstructed track in the ID are also removed. 

Events are classified in the electron (muon) channel by the presence of one electron 
(muon) with px > 25 GeV, no additional electron or muon with px > 15 GeV and at least 
four reconstructed jets with px > 25 GeV and |p| < 2.5 where at least two are identified as 
6-jets. To reduce the background contribution from non-prompt or fake leptons, E™'®® > 
30 GeV and mx(lV) > 35 GeV are also required. To reduce the effects of jet merging and 
migrations within the px ordering of the jets described in section 2, events with a pair of 
reconstructed jets separated by AR < 0.5 are vetoed. 
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5.3 Estimation of backgrounds 

The dilepton tt final states constitute the most important background to this analysis, 
followed by single top-quark production, by W boson production in association with jets 
including jets initiated by charm- and bottom-quarks iW -l-jets) and by processes where 
only many jets are produced (multi-jet production). In comparison, Z boson production 
in association with jets (Z-|-jets) and diboson production processes are smaller background 
components. 

The fraction of dilepton events that remain after applying the full event selection is 
evaluated with a MC simulation and removed from the data. A bin-by-bin correction 
factor derived from the baseline powheg-I-pythia ti Monte Carlo simulation is used (see 
section 8). Contributions from single top-quark, Z-|-jets and diboson iWW, WZ, ZZ) 
production are evaluated using corresponding MC samples and theoretical cross-sections 
for these processes, as discussed in section 4. 

The overall normalisation of the W -|-jets MC sample was determined in the data via 
a lepton charge asymmetry measurement described in ref. [15]. This method exploits the 
fact that the production of W bosons at the LHC is charge-asymmetric and also that 
the ratio of the number of W~ to bosons is more precisely known than the total 
number of W bosons [80]. Most of the other background processes result in lepton charge 
distributions that are symmetric. The numbers of events with positively and negatively 
charged leptons are measured in the data and are referred to as and respectively. 
A MC simulation was used to estimate the charge-asymmetric background from single top- 
quarks and to subtract that from the values of and The number of VL-l-jets events 
was then extracted from: 

+ Ny,- = - Ni) 

rue - I 

where tmc is the ratio of the IT "’'-l-jets and IT“-|-jets production cross-sections determined 
using the IT-|-jets MC simulation for the signal region kinematic cuts, and Nyr+ (A^^_) 
is the number of IT"*" {W~) events. The W -|-jets normalisation was determined using 
the event selection of this analysis, but without the 6-tagging requirement. The values of 
{Nyr+ + Nyr-) are independently determined for VC-|-4-jet and W+ > 5-jet events and 
are applied as inclusive scale factors to the alpgen VL+jet processes. The normalisation 
was separately obtained for each of the MC systematic uncertainty evaluations listed in 
section 6. 

The normalisation of the heavy-flavour fractions within the W -|-jets sample was deter¬ 
mined by measuring the number of IT-|-2-jet events, without a 6-tagging requirement and 
with the requirement of at least one 6-tag. To make this measurement, the charge asymme¬ 
try technique described above was applied to both sets of selected events. The number of 
events that have one or more 6-tags is related to the number of events before the 6-tagging 
requirement, the 6-tagging probability and the flavour fractions in the IT-|-jets sample. The 
simulated fractions for the heavy-flavour processes (IT66-|-jets, ITcc-|-jets and ITc-|-jets) 
are determined with the ratio WccjWbb taken from simulation. The measurement used 
the number of events after requiring at least one 6-tag where the overall normalisation of 
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VF+jets events was fixed using the values previously determined by the charge asymmetry 
method. The heavy-flavour fractions are then extrapolated from the IT-|-2-jet selection to 
W+ > 4-jet events using the heavy-flavour fractions from the MC simulation. As a result 
the Wbb+iets process is scaled up by about 30% with respect to the Monte Carlo prediction 
with a relative systematic uncertainty of 20% and the VTc-l-jets process is scaled down by 
about 15% with a relative systematic uncertainty of 40%. The scale factors are determined 
separately for the electron and muon channel. 

Multi-jet production processes have a large cross-section and can provide a non-prompt 
or fake lepton signature due to interactions with detector material, electromagnetic shower 
fluctuations and heavy-flavour decays. In the electron channel, jets and electrons from 
photon conversions or heavy-flavour decays can mimic isolated electrons from W bosons. 
In the muon channel, the background is dominated by the decay of heavy-flavour hadrons 
to muons. 

A matrix method [81] is used to estimate the number of background events using a 
second event sample for which the lepton identihcation criteria are relaxed and the isolation 
requirements removed (loose selection). The number of background events that pass the 
standard tight lepton selection described in section 5.2 is then given by: 

^Ught ^ -ffake-- ATtigl'*), 

^real ^fake 

where Creai a-iid €fake are the fractions of real and fake leptons that pass the loose and the 
tight selection and is the number of events with a lepton passing the loose 

(tight) selection. The efficiency for a real lepton to pass the tight selection (creai) is mea¬ 
sured using a tag-and-probe technique using leptons from Z boson decays. The efficiency 
for a loose background event to pass the tight selection (efake) is estimated in control re¬ 
gions dominated by background. Contributions from IT-|-jets and Z-|-jets production are 
subtracted in the control regions using simulation. 

In the electron channel the control region to determine Cfake is defined by < 

20 GeV. In the muon channel two control regions with similar numbers of events are used. 
One control region has events with low-mT(IT) while the other control region contains 
events where the selected muons have a large impact parameter. The efficiency Cfake is 
extracted separately from the two background-enriched samples and the average is used. 

6 Systematic uncertainties 

The systematic uncertainties due to detector effects and the modelling of signal and back¬ 
ground are determined for each bin of the measured observables. Each systematic un¬ 
certainty is evaluated by varying the relevant source by one standard deviation about its 
nominal value. This effect is propagated through the event selection, unfolding and correc¬ 
tion procedure. Deviations from the nominal case are evaluated separately for the upward 
and downward variations for each bin, observable and channel. The total systematic un¬ 
certainty for each bin is calculated by adding the individual systematic contributions for 
that bin in quadrature. 
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The uncertainty on the tt signal and background modelling are discussed in sections 6.1 
and 6.2. The detector-related uncertainties are presented in section 6.3. 

6.1 Signal modelling 

The uncertainty due to the choice of MC generator to model ti events is estimated by 
comparing MC@NLO+HERWlG and ALPGEN+HERWIG, while the colour reconnection mod¬ 
elling uncertainty is estimated by comparing the nominal powheg-|-pythia sample to the 
associated MC tune where the colour reconnection is disabled in pythia. 

The uncertainty in modelling the parton shower and hadronisation is evaluated from 
the relative difference between the alpgen-|-pythia tt MC sample and the alternative 
ALPGEN-I-HERWIG Sample. 

The evaluation of the uncertainty due to the choice of PDF set is obtained using the 
NNPDF 2.0 [82], mstw2008nlo and GTEQ66 [83] PDF sets. An envelope of uncertainty 
bands is determined using the PDF4LHC recommendations [84]. 

The uncertainty associated with the modelling of additional QCD radiation accompa¬ 
nying the tt system is calculated by comparing the alpgen-|-pythia sample to the ones 
with varied radiation settings presented in section 4. The variation is achieved by changing 
the renormalisation scale associated with ag consistently in the hard-scattering matrix ele¬ 
ment as well as in the parton shower. The level of radiation through parton showering [85] is 
adjusted to encompass the ATLAS measurement of additional jet activity in tt events [30]. 
The uncertainty is estimated as the maximum difference between the specialised samples 
and the nominal sample, with the uncertainty being symmetrised. 

6.2 Background modelling 

The individual experimental and theoretical uncertainties are used to calculate the uncer¬ 
tainty on the size of background contributions determined by MC simulation. This results 
in an uncertainty on the background subtraction for all backgrounds except for those from 
W -fjets processes and from multi-jet processes resulting in a non-prompt or fake lepton 
signature. 

In the muon channel the normalisation and shape uncertainties are determined as the 
difference between the two multi-jet estimates described in section 5.3. The normalisation 
uncertainty is determined to be 20%. The shape uncertainty is evaluated by two different 
linear combinations of the two multi-jet estimates. In the electron channel the normali¬ 
sation uncertainty is determined from the variation of the background estimate when the 
efficiencies for real and fake leptons are varied within their uncertainties. The uncertainty 
of the efficiency for real leptons is estimated by varying the fit parameters in the tag-and- 
probe method using Z boson events. The uncertainty on the efficiency for fake leptons is 
estimated by varying the cut between 15 GeV and 25 GeV and by relaxing the 6-tag 

requirement to at least one 6-tag. From these variations a normalisation uncertainty of 
50% is assigned. 

For the background contribution due to kF-f-jets, the overall uncertainty from the charge 
asymmetry normalisation method (see section 5.3) and the uncertainty on the flavour frac¬ 
tions are separately determined. A shape uncertainty is estimated by varying model pa- 
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rameters in the hh+jets alpgen simulation. The total background uncertainty is evaluated 
by adding in quadrature each of the different background uncertainty contributions. 

6.3 Experimental uncertainties 

The experimental uncertainties refer to the quality of the detector simulation to describe 
the detector response in data for each of the reconstructed objects. These uncertainties 
affect the MC signal and background predictions, changing the numbers of events accepted. 

The jet energy scale (JES) systematic uncertainty [77] is a major contributor to the 
overall systematic uncertainty in all distributions affected by the signal efficiency and bin 
migration. In the central region of the detector (jryj < 1.7) it varies from 2.5% to 8% 
as a function of jet pT and rj as estimated from in situ measurements of the detector 
response [77]. It incorporates uncertainties from the jet energy calibration, calorimeter re¬ 
sponse to jets, detector simulation, and the modelling of the fragmentation and underlying 
event, as well as other choices in the MC generation. Additional sources of the JES uncer¬ 
tainty are also estimated. The main contributions are the intercalibration of the forward 
region detector response from the central regions of the detector, effects from the correction 
of additional pp interactions, jet flavour composition, 6-jet JES calibration and the pres¬ 
ence of close-by jets. Uncertainties due to different detector-simulation configurations used 
in the analysis and in the calibration are added as one additional uncertainty parameter 
(“relative non-closure”). The JES uncertainty is evaluated using a total of 21 individ¬ 
ual components to model the uncertainty correlations as a function of the jet transverse 
momentum and the rapidity. 

The jet energy resolution (JER) has been found to be well modelled in simulation. In 
situ methods are used to measure the resolution, which MC simulation describes within 
10% for jets in the px range 30-500 GeV [86]. The jet reconstruction efficiency is also well 
modelled by the simulation and the uncertainty is evaluated by randomly removing simu¬ 
lated jets within the Icr uncertainty of jet reconstruction efficiency measured in data [77]. 

The uncertainties introduced by the JVF requirement used to suppress pile-up jets are 
estimated using events with a leptonic Z boson decay and an associated high-px jet. The 
efficiency to select jets from the hard-scatter and the contamination by jets produced by 
pile-up interactions are measured in appropriate control regions and the agreement between 
data and MC simulation is evaluated. 

The uncertainties related to the MC modelling of the lepton trigger, reconstruction and 
identihcation efficiency are evaluated by comparing high-purity events featuring leptons in 
data and simulation. These include Z —)• ee, Z —> pp and VU —)■ ei^ events in data 
and simulation, while tt events are also included in the simulation studies [73]. Similar 
studies are also performed for the lepton energy and momentum scales and resolutions. 
Since the two channels require different lepton flavours, the electron uncertainty affects 
only the electron channel, and similarly for the muon channel. The electron uncertainty is 
approximately double the muon uncertainty. In both cases the uncertainty is small with 
little variation between bins. The uncertainty on E™'®® is determined by propagating all 
the uncertainties associated with the energy scales and resolutions for leptons and jets to 
the calculation of E™'®®. Two additional uncertainties that are included originate from 
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the calorimeter cells not associated with any physics objects and the pile-up modelling. 
The systematic uncertainties associated with tagging jets originating from 6-quarks are 
separated into three categories. These are the efficiency of the tagging algorithm (6-quark 
tagging efficiency), the efficiency with which jets originating from c-quarks pass the 6- 
tag requirement (c-quark tagging efficiency) and the rate at which light-flavour jets are 
tagged (misidentified tagging efficiency). The efficiencies are estimated from data and 
parameterised as a function of px and rj [78, 87]. The systematic uncertainties arise from 
factors used to correct the differences between simulation and data in each of the categories. 
The uncertainties in the simulation modelling of the 6-tagging performance are assessed 
by studying 6-jets in dileptonic ti events [88]. The 6-tagging efficiency is another major 
contributor to the overall systematic uncertainty and tends to slightly increase with px- The 
c-quark efficiency uncertainty is approximately constant at ~ 2% while the misidentification 
uncertainty contributes at the percent level for all distributions. 

The last experimental uncertainty evaluated is due to the measurement of the inte¬ 
grated luminosity. This is dominated by the accuracy of the beam separation scans and 
has an associated uncertainty of 1.8% [72] that is assigned to each bin of the distributions 
and the MC background predictions. With the exception of the MC and data statistics 
and the background modelling all uncertainty components are correlated across the bins 
and for all observables. 

7 Reconstructed yields and distributions 

A summary of the number of selected data events, background contributions and total 
predictions is given in table 1. Dilepton tt events constitute the largest background followed 
by single top-quark production. The lT-|-jets and non-prompt or fake lepton backgrounds 
are smaller in comparison. 

The reconstructed distributions of the observables of the t\, th and t\t\i system are shown 
in figures 3 to 7. This includes the px and rapidity of the individual pseudo-top-quarks 
and the px, rapidity and mass of the tith system for the muon and electron channels. 

In each region, the expected number of events agrees with the number observed in 
the data. The data are shown using bin sizes that correspond to one standard deviation 
resolution in the hadronic variables, except in the tails of the distributions where the bin 
width is increased to reduce the statistical fluctuations. 
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Yields 

Source 

e -b jets (ibstat. ± syst.) 

± jets (±stat. ± syst.) 

tt (single-lepton) 

7810 ± 20 ± 630 

9410 ± 22 ± 700 

tt (dilepton) 

570 ± 5± 60 

670 ± 6± 60 

Single top-quarks 

430 ± 6± 50 

520 ± 7± 50 

Vk-bjets 

260 ±15 ±110 

360 ± 18 ± 130 

Non-prompt or fake leptons 

130±16± 70 

56 ± 4± 11 

Z-bjets 

52 ± 4± 13 

28 ± 3 ± 7 

Diboson 

7± 1 ± 1 

7± 1± 1 

Expectation 

9270 ± 31 ± 650 

11050 ±30 ±720 

Observed 

8791 

10690 


Table 1. The number of selected events from data, MC simulation and data-driven background 
estimates, in the electron and muon channels. The yields of MC samples that are not constrained by 
fits to data shown are normalised to an integrated luminosity of 4.6 fb“^. The uncertainties specific 
to the ti MC sample (ISR, generator, fragmentation, etc.) are not included in the tt (single-lepton 
and dilepton) uncertainties given in this table. 
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Figure 3. The reconstructed pseudo-top-quark px in comparison to the MC signal and data-driven 
background models. The PT(th) distributions: (a) the muon channel and (b) the electron channel. 
The PT(ti) distributions: (c) the muon channel and (d) the electron channel. Signal and background 
processes are shown in different colours, with “Other” including the small backgrounds from diboson 
and Z+]ets production, and non-prompt or fake lepton signatures from multi-jet processes. The 
data are compared with predictions from background models and expected yields from simulated 
ti events generated using powheg-|-pythia with the “C” variant of the Perugia 2011 tunes family. 
The error bars on the data points show the statistical uncertainty on the data, while the shaded 
band shows the total systematic and statistical uncertainty on the predicted yields. The systematic 
uncertainty has a strong bin-to-bin correlation. 
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Figure 4. The reconstructed pseudo-top-quark rapidity in comparison to the MC signal and 
data-driven background models. The |y(ih)| distributions: (a) the muon channel and (b) the elec¬ 
tron channel. The \y{t\)\ distributions: (c) the muon channel and (d) the electron channel. The 
POWHEG-I-PYTHIA MC generator with the Perugia 2011C tune is used for the tt signal estimate. 
Signal and background processes are shown in different colours, with “Other” including small back¬ 
grounds from diboson and Z-f-jets production, as well as non-prompt and fake leptons from multi-jet 
processes. The shaded band shows the total systematic and statistical uncertainties on the signal 
plus background expectation. 
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Figure 5. The reconstructed pt for the system of pseudo-top-quark pairs iii^ in comparison to the 
MC signal and data-driven background models. The PT{tith) distributions: (a) the muon channel 
and (b) the electron channel. Signal and background processes are shown in different colours, with 
“Other” including small backgrounds from diboson and Z-|-jets production, as well as non-prompt 
and fake leptons from multi-jet processes. The powheg-|-pythia MC generator with the Perugia 
2011C tune is used for the tt signal estimate. The shaded band shows the total systematic and 
statistical uncertainties on the signal plus background expectation. 
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Figure 6. The reconstructed rapidity for the pseudo-top-quark system iii^ in comparison to the 
MC signal and data-driven background models. The y(tith) distributions: (a) the muon channel 
and (b) the electron channel. Signal and background processes are shown in different colours, with 
“Other” including small backgrounds from diboson and Z-|-jets production, as well as non-prompt 
and fake leptons from multi-jet processes. The powheg-|-pythia MC generator with the Perugia 
2011C tune is used for the tt signal estimate. The shaded band shows the total systematic and 
statistical uncertainties on the signal plus background expectation. 


- 23 - 




















































(a) (b) 

Figure 7. The reconstructed TO(tith) for the pseudo-top-quark system tith in the muon (a) and 
electron (b) channel. The distributions are shown in comparison to the MC signal and data-driven 
background models. Signal and background processes are shown in different colours, with “Other” 
including small backgrounds from diboson and Z-|-jets production, as well as non-prompt and fake 
leptons from multi-jet processes. The POWHEG-|-pythia MC generator with the Perugia 2011C 
tune is used for the tt signal estimate. The shaded band shows the total systematic and statistical 
uncertainties on the signal plus background expectation. 
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8 Corrections for detector effects and channel combinations 

Each of the reconstructed pseudo-top-quark observable distributions is corrected for the 
effects of detector efficiencies and resolution. All distributions are presented within the 
kinematic range dehned in section 2, which is close to the acceptance of the reconstructed 
object and event selections, such that model dependencies from regions of phase space 
outside of the acceptance are minimised. Section 8.1 describes the correction procedure, 
section 8.2 describes the propagation of the statistical and systematic uncertainties to the 
final distributions, and section 8.3 describes the combination of the results obtained in the 
electron and muon channels. 


8.1 Correction procedure 

The reconstructed pseudo-top-quark observable distributions are corrected as follows: 


AT* = f* , • . P . 

-‘'part ./partlreco recoj J misassign •trecolpart V'reco 

3 



where Nieco (-^part) is the number of reconstructed (fully corrected) events in a given 
reconstructed observable bin j (particle-level observable bin i), and is the number of 

background events estimated as explained in section 5.3. The correction factors /partlreco’ 
/misaasign ^nd /.^ecolpart ^re disCUSSed below. 

Detector resolution effects on the reconstructed pseudo-top-quark observables are cor¬ 
rected with an iterative Bayesian unfolding procedure [89] using a response matrix Mpa™ 
that describes the migration between the detector-level observable bin j and particle- 
level observable bin i. To ensure a one-to-one relationship between the particle-level and 
the detector-level observables, the response matrix is constructed from events where for 
each detector-level pseudo-top-quark object a particle-level pseudo-top-quark object can 
be matched. The matching is based on the angular differences between the components 
of the pseudo-top-quarks. Detector-level and particle-level jets (jet/jet') and leptons (i/i') 
are considered matched if they satisfy Ai2(jet, jet') < 0.35 and < 0.02. This 

matching definition is found to be fully efficient for leptons and close to 100% efficient for 
jets. 

The response matrix is derived from MC simulations and has diagonal elements that 
are larger than 0.6 for all observables. The elements in the matrix are normalised to the 
total number of reconstructed events determined from the sum of all the bins associated 
with the x-axis. The response matrix is applied using two iterations, which is found to 
provide convergence and avoid higher statistical uncertainties in the tails of the corrected 
distributions. 

Events that have no matched detector-level and particle-level pair are taken into ac¬ 
count by three factors that correct the detector-level observable distribution to the particle- 
level observable distribution: 


• The correction for events that pass the detector-level event selection but fail the 
particle-level event selection (/recolpart); 
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• The correction factor for events with a reconstructed pseudo-top-quark that has no 
counterpart at the particle level (/misassign); 

• The correction for events that fulfil the particle-level event selection requirements but 
fail the reconstruction-level event selection (/partireco)- 

These correction factors are also derived from MC simulation. The correction factors 
/recolpart and /misassign are in the range between 0.65-0.70 for all observables. The cor¬ 
rection factor /partireco is primarily dominated by the detector efficiency, in particular by 
the 6-tagging efficiency, and is in the range 5.5-8.0. They are found to be similar for the 
ti generators discussed in section 4, such that no large MC modelling dependencies are 
observed between the different MC samples. 

The same unfolding procedure is applied to each of the distributions PT{t), y{t), y{tith), 
PT{iiih), and m{iiih). 

The number of background events A^bgnd and the correction factors /recolpart and 
/misassign are functions of the detector-level pseudo-top-quark observables. The correction 
factor /partireco is a function of the particle-level pseudo-top-quark observable Xp^rt- To 
evaluate the cross-section in bin i, it is also necessary to take account both the luminosity 
and bin width. 

8.2 Propagation of uncertainties 

With the exception of the non-tt backgrounds, each of the correction factors in eq. (8.1) is 
calculated using powheg-|-pythia events that are passed through the detector simulation. 
The effect of the statistical uncertainty on Mp“° is estimated by smearing the number of 
events in each element of the matrix, using a Poisson probability density function. The sta¬ 
tistical uncertainty on the correction factors (/recolpart) /misassign and /partireco) is evaluated 
by smearing the value in each bin using a Gaussian distribution following the statistical 
uncertainty in the bin. The correction factors and the response matrix are smeared simul¬ 
taneously by performing 1000 pseudo-experiments and repeating the unfolding procedure 
for each pseudo-experiment. The statistical uncertainty for each measurement point is 
taken from the root-mean-square (la) of the spread of the unfolded distributions over the 
various pseudo-experiments. 

The statistical uncertainty on the reconstructed distributions in data is propagated 
to the hnal distributions by performing 1000 pseudo-experiments, following a Poisson dis¬ 
tribution defined by the number of events in each bin j. Similar to the MC statistical 
uncertainty, each bin of xleco is independently fluctuated. 

The experimental systematic uncertainty on the reconstructed distributions is evalu¬ 
ated by changing the values of the physics objects by their associated uncertainties. The 
total uncertainty on the number of reconstructed background events (Wgnd) is evaluated 
by summing in quadrature each of the background uncertainties discussed in section 6. 

The systematic uncertainty on the unfolded spectra due to the background is evalu¬ 
ated by performing 1000 pseudo-experiments, following a normal distribution with a width 
matching the total uncertainty band. The root-mean-square of the distribution of un¬ 
folded spectra of the pseudo-experiments is taken as the uncertainty on the background. 
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Due to the small number of background events and their generally small uncertainties, the 
background uncertainties are propagated in an uncorrelated manner. This choice reduces 
the impact of statistical fluctuations on the background uncertainty estimates for the final 
distributions. 

The systematic uncertainties on /recolpart, /misassign and /partlreco arising due 

to the choice of the powheg+pythia ti MC model are each evaluated as a relative bias. 
For a given MC simulation set the bias is dehned as the fully corrected unfolded yield 
for a given bin minus the true particle yield for that bin. The relative bias is defined as 
the difference in bias between the nominal MC sample and the varied MC sample, both 
unfolded using the nominal MC sample. For the powheg+pythia ti MC sample, the bias 
is found to be consistent with zero within the statistical uncertainties of each measurement 
point. For each ti modelling systematic uncertainty, a pair of particle-level and detector- 
level spectra is generated. One thousand pseudo-experiments are used to fluctuate the 
reconstructed input spectrum within its statistical uncertainty. The pseudo-experiments 
are used to evaluate the statistical significance of the systematic variation in the output 
distribution. The relative bias is calculated for each pseudo-experiment. 

The ISR/FSR systematic uncertainty is evaluated from the relative bias between the 
ALPGEN+PYTHIA central and shifted ISR/FSR samples. The uncertainty on the matrix- 
element calculation and matching scheme (the generator uncertainty) is estimated from the 
relative bias of mc@nlo+herwig with respect to the alpgen+herwig ti sample. 

Each of the ti uncertainties is propagated individually and are then symmetrised before 
being combined, taking the larger of the upward or downward variation. 

8.3 Combination of lepton channels 

The electron and muon channel measurements of each pseudo-top-quark distribution are 
combined by using the Best Linear Unbiased Estimate (BLUE) method [90, 91]. The 
BLUE method determines the coefficients (weights) to be used in a linear combination of 
the input measurements by minimising the total uncertainty of the combined result. All 
uncertainties are assumed to be distributed according to a Gaussian probability density 
function. The algorithm takes both the statistical and systematic uncertainties and their 
correlations into account. The BLUE combination was cross-checked against an average 
performed using the algorithm discussed in ref. [92]. The result of the two methods are 
found to be consistent. The MC statistical uncertainties on the correction factors for the 
two samples are assumed to be uncorrelated. The uncertainties related to the electron 
and muon efficiencies are also treated uncorrelated. All other systematic uncertainties 
are treated as fully correlated. In particular, the total background systematic uncertainty 
is assumed to be completely correlated between the electron and muon channel. The 
uncertainties of the combined measurement are dependent on the observables but tend to 
closely follow the muon channel uncertainties. 


- 27 - 


9 Results 


The measurements of the differential tt cross-section corrected for detector effects are pre¬ 
sented as a function of the PT{t) and y{t) for the hadronic or leptonic reconstruction of 
the pseudo-top-quarks, as well as the variables yitit^), PT{tith) and m(tith) of the pseudo- 
top-quark-pair system; t refers to both the hadronic and leptonic pseudo-top-quark. The 
fiducial cross-section measurements are presented within the kinematic range defined in 
section 2, and are evaluated as described in section 8. 

The measurements can be directly compared to Monte Carlo simulations that are 
generated with top quarks with a mass of 172.5 GeV. The top quark mass dependence of 
the cross-sections was derived with Monte Carlo samples generated with top quark masses 
from 165 to 178 GeV. The cross-sections depend linearily on the top mass. The slope® is 
between —0.1 and —0.3 %/GeV for y{t), y(tith) and PT{tith)- For m(tith) (pt(c) it rises 
from —0.6 (—0.3) %/GeV for low values to 0.4 (0.4) %/GeV at high values. 

The tt MC generators are of two types (for a detailed discussion see section 4): NLO 
matrix elements are used to describe the hard scattering pp —)• ft, and LO multi-leg MC 
generators model ti processes with up to five additional quark or gluon emissions. All ti 
MC samples are normalised to the NNLO-I-NNLL QCD inclusive cross-section. 

Figures 8, 9 and 10 show the corrected data for the differential variables noted above. 
Superimposed on the data are the expectations of the LO multi-leg MC generator alpgen 
(see section 4): 

• ALPGEN interfaced with pythia (alpgen-|-pythia), 

• ALPGEN interfaced with herwig (alpgen-|-herwig), 

• ALPGEN interfaced with pythia but with the renormalisation scale varied by a factor 
of 2 (as Up) and a factor of 0.5 (as Down). 

Both the ALPGEN-I-HERWIG and alpgen-|-pythia models indicate yields within the 
acceptance that are higher than the observed data. This is evident in figures 9 and 10(b) 
where for the complete rapidity range the models predict more events. The models also 
indicate a PT{t) spectrum that is harder than the ones observed in data, a possible con¬ 
sequence of using the GTEQ6l1 PDF set. For the alpgen-|-pythia sample, the effect of 
increased or decreased radiation is illustrated with the renormalisation scale changed by 
a factor of 2.0 (as Up) and by a factor of 0.5 (as Down) applied consistently to both 
ALPGEN and PYTHIA as noted in section 4. The increased radiation gives fewer events at 
low pt(d and more at high pt(C- This is at the level of 5-10%. For the leptonic PT{t) 
this effect is slightly larger as shown in figure 8(b). The PT(b^h) distribution, shown in 
figure 10(a), is very sensitive to additional radiation. When the renormalisation scale fac¬ 
tor changes from 0.5 to 2.0, 15-20% more events are observed at low PT(b^h) and 20-30% 
fewer at high PTityth)- Increased radiation also leads to 5% fewer events at low m(tith) and 
10% more events at high m(tith) as shown in figure 10(c). Nevertheless, the (as Down) 

®The top mass (rut) dependence of the cross-sections can be obtained from the slope a quoted above via 
a-{mt- 172.5). 


- 28 - 



CL 


10 ' 


M 0'“ 


10'^ 


10-' 


ATLAS is = 7 TeV, 4.6 fb ’ 

—•— Data 

' „ stat. 0 syst. uncart. 

™ ALPGEN+PYTHIA 

* - ALPGEN+PYTHIA (a^ UP) 

Hp . ALPGEN+PYTHIA (ct^ DOWN) 

^ ALPGEN+HERWIG 



CL 


10' 


0 - 10 '“ 

T3 


10'^ 


ATLAS is = 7 TeV, 4.6 fb ’ 

—•— Data 

stat. 0 syst. uncart. 

^ ALPGEN+PYTHIA 

^ ^ - ALPGEN+PYTHIA (a^ UP) 

*** . ALPGEN+PYTHIA (a^ DOWN) 

.M.M.. ALPGEN+HERWIG 


Q 

O 

o 

o 

Q. 

X 

LU 


BfFirmrrFTTmTrrFKrFTTmrrFTTOTTrFi 


0-^ 


1.8 


1.6 


1.4 


1.2 


1 


0.8 


0.^ 



200 400 600 800 

Pt(!,) [GeV] 


(a) 


(b) 


Figure 8. Differential ti cross-section after channel combination as a function of (a) the hadronic 
pseudo-top-quark pT(ih) and (b) the leptonic pseudo-top-quark pt(^i)- The data points are shown 
with a blue band which represents the total uncertainty (statistical and systematic). The model 
predictions from several LO multi-leg MC generators described in the text are superimposed: the 
ALPGEN-I-HERWIG and the ALPGEN-I-PYTHIA samples. For the ALPGEN-|-pythia sample, variations 
in the renormalisation scale by factors of 0.5 and 2.0 are also shown. 


variation is not sufficient to restore agreement between the data and the MC simulation. 
The ALPGEN-I-HERWIG Sample follows the alpgen-|-pythia sample with the (as Down) 
variation. 

Figures 11, 12 and 13 compare the data with expectations of the NLO MC generators 
POWHEG and MG@NLO. In particular, the following NLO variants are shown: 

• POWHEG- hvq V.4 generator with the gtIOnlo PDF set interfaced with pythia using 
the “C” variant of the Perugia 2011 tunes and GTEQOlI PDF set (powheg-|-pythia). 

• POWHEG- hvq v.4 generator where the gtIOnlo PDF set is replaced with the 
herapdf15nlo pdf set to assess the sensitivity of the distributions to changes in 
the gluon PDF (powheg(herapdf)-I-pythia). 

• MG@NLO generator with the gtIOnlo PDF set interfaced with herwig and jimmy 
with the AUET2 tune (mg@nlo-|-herwig). 

The individual hadronic and leptonic y{t) distributions in figure 12 are well described 
by each of the NLO MC models. The powheg(herapdf)-|-pythia sample is the one 
closest to predicting the data, whereas the other NLO models predict a more forward 
distribution. This is also true for the hadronic and leptonic PT{t) distributions (see figure 
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Figure 9. Differential ti cross-section after channel combination as a function of (a) the hadronic 
pseudo-top-quark rapidity y(th) and (b) the leptonic pseudo-top-quark rapidity y(ti). The data 
points are shown with a blue band which represents the total uncertainty (statistical and system¬ 
atic). The model predictions from several LO multi-leg MC generators described in the text are 
superimposed: the ALPGEN-|-herwig and the ALPGEN-|-pythia samples. For the ALPGEN-|-pythia 
sample, variations in the renormalisation scale by factors of 0.5 and 2.0 are also shown. 


11) for which the powheg(herapdf)-|-pythia sample lowers the cross section at high 
PT{t) with respect to the nominal sample, which results in a better description of the 
data. For these variables MC@NLO also gives a good description. ThepT(iiih) distribution, 
shown in figure 13(a), highlights the different hard-gluon emission models. The mc@nlo 
prediction is lower than the data at high PTitith), as expected due to the softer fifth-jet px 
from MC@NLO in comparison to the other generators [31]. 

All the models predict an excess in the lower m(tith) region, as shown in figure 13(c), 
implying that the threshold region description is inadequate. The powheg(herapdf)-|-pythia 
sample agrees well with the high-mass tail, while the other samples overestimate 

the tail. This is consistent with the softer gluon component in the herapdf15nlo PDF set 
compared to the one in the ctIOnlo PDF set. The 2/(tith) distribution is reasonably pre¬ 
dicted for all models in the low y(tith) region, but only the powheg-|-pythia model with 
the herapdf15nlo PDF set provides a good overall description. ThepT(iiih) spectrum is 
sensitive to the extra radiation produced in the parton collision process. All models agree 
with the data within the systematic uncertainties. If these uncertainties can be reduced, 
this suggests that the PT{tith) distribution can be used to constrain phenomenological ra¬ 
diation parameters in future MC tunes. 
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Figure 10. Differential tt cross-section after channel combination as a function of (a) the total 
leptonic and hadronic tt pseudo-top-quark variable PT(iffh), (b) the rapidity y(iiih) and (c) the 
mass m(tith). The data points are shown with a blue band which represents the total uncertainty 
(statistical and systematic). The model predictions from several LO multi-leg MC generators de¬ 
scribed in the text are superimposed: the ALPGEN-|-herwig and the ALPGEN-|-pythia samples. For 
the ALPGEN-bPYTHiA sample, variations in the renormalisation scale by factors of 0.5 and 2.0 are 
also shown. 
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Figure 11. Differential tt cross-section after channel combination as a function of (a) the 
hadronic pseudo-top-quark pT(ih) and (b) the leptonic pseudo-top-quark The data 

points are shown with a blue band which represents the total uncertainty (statistical and sys¬ 
tematic). The model predictions from several NLO MC generators described in the text are 
superimposed: powheg(ct10)-I-pythia, powheg(herapdf)-|-pythia, powheg-I-herwig and 
MC@NLO-|-HERWIG. 
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Figure 12. Differential tt cross-section after channel combination as a function of (a) the hadronic 
pseudo-top-quark rapidity y{th) and (b) the leptonic pseudo-top-quark rapidity y{ti). The data 
points are shown with a blue band which represents the total uncertainty (statistical and sys¬ 
tematic). The model predictions from several NLO MC generators described in the text are 
superimposed: powheg(ct10)-I-pythia, powheg(herapdf)-|-pythia, powheg-I-herwig and 
MC@NLO-|-HERWIG. 
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Figure 13. Differential ti cross-section after channel combination as a function of (a) the to¬ 
tal leptonic and hadronic ti pseudo-top-quark variables PT(idh), (b) the rapidity y{t\t\i) and (c) 
the mass The data points are shown with a blue band which represents the total un¬ 

certainty (statistical and systematic). The model predictions from several NLO MC generators 
described in the text are superimposed: powheg(ct10)-|-pythia, powheg(herapdf)-|-pythia, 
POWHEG-I-HERWIG and MC@NLO-|-HERWIG. 
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10 Conclusions 


Differential fiducial tt cross-section measurements are presented for kinematic variables of 
the pseudo-top-quark (t), defined at the particle level by the decay products of the W 
boson and 6-quark occurring in top-quark decays. 

The pseudo-top-quark approach is a new tool to probe QCD in the top-quark sector. 

It is an experimental observable that is correlated with the top-quark parton and is used to 
define differential tt cross-sections with reduced model dependence. It can also be used to 
assess how well MC simulations can describe the tt production mechanism in proton-proton 
collisions and to compare various MC models. 

The present measurements were performed within a kinematic range that closely 
matches the acceptance of the reconstructed objects and each reconstructed kinematic 
variable distribution is corrected (unfolded) for the effects of detector efficiency and reso¬ 
lution. 

The differential fiducial cross-sections are measured with the ATLAS detector at the 
LHC for proton-proton collisions at = 7 TeV and an integrated luminosity of 4.6 fb“^, as 
a function ofpT(ih )5 y{th), PTiti), vitith), PT{tith), and m{tith), where ti and th refers 
to the hadronic or leptonic pseudo-top-quark. The distributions provide complementary in¬ 
formation and show some sensitivity to the selected PDF set, the parton shower procedure, 
and to a lesser extent the matrix element and parton shower matching scheme. The larger 
acceptance for the alpgen MC models resulted in an excess of events reconstructed within 
the fiducial region. The higher as variation is seen to be disfavoured at higher PT(iffh) and 
m(tith) values. Among the several tt MC models used, the powheg(herapdf)-|-pythia 
sample provides the best representation for all of the distributions, except for low m{tith) 
where many of the MC models predict a higher cross-section than the data at low mass 
values. The mc@nlo prediction is seen to produce a recoil distribution that is too soft with 
respect to the data. These measurements are currently limited by the systematic uncer¬ 
tainty, the main components being the 6-tagging uncertainty, the jet energy measurement 
uncertainty and the modelling uncertainty of the initial and final state parton showers. 
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Institut fiir Experimentelle Physik IV, Technische Universitat Dortmund, Dortmund, 
Germany 

Institut fiir Kern- und Teilchenphysik, Technische Universitat Dresden, Dresden, 
Germany 

Department of Physics, Duke University, Durham NG, United States of America 
SUPA - School of Physics and Astronomy, University of Edinburgh, Edinburgh, United 
Kingdom 

INFN Laboratori Nazionali di Frascati, Frascati, Italy 

Fakultat fiir Mathematik und Physik, Albert-Ludwigs-Universitat, Freiburg, Germany 
Section de Physique, Universite de Geneve, Geneva, Switzerland 

50 (a) Sezione di Genova; Dipartimento di Fisica, Universita di Genova, Genova, 

Italy 

51 (“) Andronikashvili Institute of Physics, Iv. Javakhishvili Tbilisi State University, 
Tbilisi; High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia 

H Physikalisches Institut, Justus-Liebig-Universitat Giessen, Giessen, Germany 
SUPA - School of Physics and Astronomy, University of Glasgow, Glasgow, United 
Kingdom 

H Physikalisches Institut, Georg-August-Universitat, Gottingen, Germany 
Laboratoire de Physique Subatomique et de Cosmologie, Universite Grenoble-Alpes, 
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CNRS/IN2P3, Grenoble, France 

Department of Physics, Hampton University, Hampton VA, United States of America 
Laboratory for Particle Physics and Cosmology, Harvard University, Cambridge MA, 
United States of America 

58 (a) Kirchhoff-Institut fiir Physik, Ruprecht-Karls-Universitat Heidelberg, Heidelberg; 

(^) Physikalisches Institut, Ruprecht-Karls-Universitat Heidelberg, Heidelberg; ZITI 
Institut fiir technische Informatik, Ruprecht-Karls-Universitat Heidelberg, Mannheim, 
Germany 

Faculty of Applied Information Science, Hiroshima Institute of Technology, Hiroshima, 
Japan 

60 (a) Department of Physics, The Chinese University of Hong Kong, Shatin, N.T., Hong 
Kong; Department of Physics, The University of Hong Kong, Hong Kong; 

Department of Physics, The Hong Kong University of Science and Technology, Clear 
Water Bay, Kowloon, Hong Kong, China 

Department of Physics, Indiana University, Bloomington IN, United States of America 
Institut fiir Astro- und Teilchenphysik, Leopold-Franzens-Universitat, Innsbruck, 
Austria 

University of Iowa, Iowa City lA, United States of America 

Department of Physics and Astronomy, Iowa State University, Ames lA, United States 
of America 

Joint Institute for Nuclear Research, JINR Dubna, Dubna, Russia 
KEK, High Energy Accelerator Research Organization, Tsukuba, Japan 
Graduate School of Science, Kobe University, Kobe, Japan 
Faculty of Science, Kyoto University, Kyoto, Japan 
Kyoto University of Education, Kyoto, Japan 
™ Department of Physics, Kyushu University, Fukuoka, Japan 

Instituto de Fisica La Plata, Universidad Nacional de La Plata and CONICET, La 
Plata, Argentina 

Physics Department, Lancaster University, Lancaster, United Kingdom 
73 (a) Sezione di Lecce; Dipartimento di Matematica e Fisica, Universita del 

Salento, Lecce, Italy 

Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 
Department of Physics, Jozef Stefan Institute and University of Ljubljana, Ljubljana, 
Slovenia 

School of Physics and Astronomy, Queen Mary University of London, London, United 
Kingdom 

Department of Physics, Royal Holloway University of London, Surrey, United Kingdom 
Department of Physics and Astronomy, University College London, London, United 
Kingdom 

Louisiana Tech University, Ruston LA, United States of America 
Laboratoire de Physique Nucleaire et de Hautes Energies, UPMC and Universite 
Paris-Diderot and CNRS/IN2P3, Paris, France 

Fysiska institutionen, Lunds universitet, Lund, Sweden 
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Departamento de Fisica Teorica C-15, Universidad Autonoma de Madrid, Madrid, 
Spain 

Institut fiir Physik, Universitat Mainz, Mainz, Germany 

School of Physics and Astronomy, University of Manchester, Manchester, United 
Kingdom 

CPPM, Aix-Marseille Universite and CNRS/IN2P3, Marseille, France 
Department of Physics, University of Massachusetts, Amherst MA, United States of 
America 

Department of Physics, McGill University, Montreal QG, Ganada 
School of Physics, University of Melbourne, Victoria, Australia 

Department of Physics, The University of Michigan, Ann Arbor MI, United States of 
America 

Department of Physics and Astronomy, Michigan State University, East Lansing MI, 
United States of America 

91 (“) INFN Sezione di Milano; Dipartimento di Fisica, Universita di Milano, Milano, 
Italy 

B.I. Stepanov Institute of Physics, National Academy of Sciences of Belarus, Minsk, 
Republic of Belarus 

National Scientific and Educational Gentre for Particle and High Energy Physics, 
Minsk, Republic of Belarus 

Department of Physics, Massachusetts Institute of Technology, Cambridge MA, United 
States of America 

Group of Particle Physics, University of Montreal, Montreal QC, Canada 
P.N. Lebedev Institute of Physics, Academy of Sciences, Moscow, Russia 
Institute for Theoretical and Experimental Physics (ITEP), Moscow, Russia 
National Research Nuclear University MEPhI, Moscow, Russia 
D.V. Skobeltsyn Institute of Nuclear Physics, M.V. Lomonosov Moscow State 
University, Moscow, Russia 

Fakultat fiir Physik, Ludwig-Maximilians-Universitat Miinchen, Miinchen, Germany 
Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut), Miinchen, Germany 
Nagasaki Institute of Applied Science, Nagasaki, Japan 

Graduate School of Science and Kobayashi-Maskawa Institute, Nagoya University, 
Nagoya, Japan 

104 (a) Sezione di Napoli; Dipartimento di Fisica, Universita di Napoli, Napoli, 

Italy 

Department of Physics and Astronomy, University of New Mexico, Albuquerque NM, 
United States of America 

106 Institute for Mathematics, Astrophysics and Particle Physics, Radboud University 
Nijmegen/Nikhef, Nijmegen, Netherlands 

Nikhef National Institute for Subatomic Physics and University of Amsterdam, 
Amsterdam, Netherlands 

Department of Physics, Northern Illinois University, DeKalb IL, United States of 
America 
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Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, Russia 

Department of Physics, New York University, New York NY, United States of America 

Ohio State University, Columbus OH, United States of America 

Faculty of Science, Okayama University, Okayama, Japan 

Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma, 
Norman OK, United States of America 

Department of Physics, Oklahoma State University, Stillwater OK, United States of 
America 

Palacky University, RCPTM, Olomouc, Czech Republic 

Center for High Energy Physics, University of Oregon, Eugene OR, United States of 
America 

LAL, Universite Paris-Sud and CNRS/IN2P3, Orsay, Prance 
Graduate School of Science, Osaka University, Osaka, Japan 
Department of Physics, University of Oslo, Oslo, Norway 
Department of Physics, Oxford University, Oxford, United Kingdom 
121 (a) jjsjpN Sezione di Pavia; Dipartimento di Eisica, Universita di Pavia, Pavia, Italy 
Department of Physics, University of Pennsylvania, Philadelphia PA, United States of 
America 

Petersburg Nuclear Physics Institute, Gatchina, Russia 
124 (a) Sezione di Pisa; Dipartimento di Eisica E. Permi, Universita di Pisa, Pisa, 

Italy 

Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh PA, 
United States of America 

126 (a) Laboratorio de Instrumentacao e Eisica Experimental de Particulas - LIP, Lisboa; 

Eaculdade de Ciencias, Universidade de Lisboa, Lisboa; Department of Physics, 
University of Coimbra, Coimbra; Centro de Eisica Nuclear da Universidade de Lisboa, 
Lisboa; Departamento de Eisica, Universidade do Minho, Braga; Departamento de 
Eisica Teorica y del Cosmos and CAFPE, Universidad de Granada, Granada (Spain); 

Dep Eisica and CEFITEC of Eaculdade de Ciencias e Tecnologia, Universidade Nova de 
Lisboa, Caparica, Portugal 

Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech 
Republic 

Czech Technical University in Prague, Praha, Czech Republic 
Faculty of Mathematics and Physics, Charles University in Prague, Praha, Czech 
Republic 

State Research Center Institute for High Energy Physics, Protvino, Russia 
131 Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United 
Kingdom 

Ritsumeikan University, Kusatsu, Shiga, Japan 

133 (a) Sezione di Roma; Dipartimento di Eisica, Sapienza Universita di Roma, 

Roma, Italy 

134 (a) Sezione di Roma Tor Vergata; Dipartimento di Eisica, Universita di 

Roma Tor Vergata, Roma, Italy 
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135 (a) Sezione di Roma Tre; Dipartimento di Matematica e Fisica, Universita 

Roma Tre, Roma, Raly 

136 (a) Paculte des Sciences Ain Chock, Reseau Universitaire de Physique des Hautes 
Energies - Universite Hassan II, Casablanca; Centre National de I’Energie des Sciences 
Techniques Nucleaires, Rabat; Faculte des Sciences Semlalia, Universite Cadi Ayyad, 
LPHEA-Marrakech; Faculte des Sciences, Universite Mohamed Premier and LPTPM, 
Oujda; Faculte des sciences, Universite Mohammed V-Agdal, Rabat, Morocco 

DSM/IRFU (Institut de Recherches sur les Lois Fondamentales de I’Univers), CEA 
Saclay (Commissariat a I’Energie Atomique et aux Energies Alternatives), Gif-sur-Yvette, 
France 

Santa Cruz Institute for Particle Physics, University of California Santa Cruz, Santa 
Cruz CA, United States of America 

Department of Physics, University of Washington, Seattle WA, United States of 
America 

Department of Physics and Astronomy, University of Sheffield, Sheffield, United 
Kingdom 

Department of Physics, Shinshu University, Nagano, Japan 

Fachbereich Physik, Universitat Siegen, Siegen, Germany 

Department of Physics, Simon Fraser University, Burnaby BC, Canada 

SLAC National Accelerator Laboratory, Stanford CA, United States of America 

145 (a) Faculty of Mathematics, Physics Sz Informatics, Comenius University, Bratislava; 
Department of Subnuclear Physics, Institute of Experimental Physics of the Slovak 

Academy of Sciences, Kosice, Slovak Republic 

146 (a) Department of Physics, University of Cape Town, Cape Town; Department of 
Physics, University of Johannesburg, Johannesburg; School of Physics, University of 
the Witwatersrand, Johannesburg, South Africa 

147 (a) Department of Physics, Stockholm University; The Oskar Klein Centre, 
Stockholm, Sweden 

148 Physics Department, Royal Institute of Technology, Stockholm, Sweden 
Departments of Physics & Astronomy and Chemistry, Stony Brook University, Stony 

Brook NY, United States of America 

Department of Physics and Astronomy, University of Sussex, Brighton, United 
Kingdom 

School of Physics, University of Sydney, Sydney, Australia 
Institute of Physics, Academia Sinica, Taipei, Taiwan 

Department of Physics, Technion: Israel Institute of Technology, Haifa, Israel 
Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 
Tel Aviv, Israel 

Department of Physics, Aristotle University of Thessaloniki, Thessaloniki, Greece 
International Center for Elementary Particle Physics and Department of Physics, The 
University of Tokyo, Tokyo, Japan 

Graduate School of Science and Technology, Tokyo Metropolitan University, Tokyo, 
Japan 
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Department of Physics, Tokyo Institute of Technology, Tokyo, Japan 
Department of Physics, University of Toronto, Toronto ON, Canada 
160 (a) tRIUMF, Vancouver BC; Department of Physics and Astronomy, York 
University, Toronto ON, Canada 

Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Japan 
Department of Physics and Astronomy, Tufts University, Medford MA, United States 
of America 

Centro de Investigaciones, Universidad Antonio Narino, Bogota, Colombia 
Department of Physics and Astronomy, University of California Irvine, Irvine CA, 
United States of America 

165 (a) j]sjp]sj Gruppo Collegato di Udine, Sezione di Trieste, Udine; ICTP, Trieste; 
Dipartimento di Chimica, Fisica e Ambiente, Universita di Udine, Udine, Italy 

Department of Physics, University of Illinois, Urbana IL, United States of America 
Department of Physics and Astronomy, University of Uppsala, Uppsala, Sweden 
Institute de Fisica Corpuscular (IFIC) and Departamento de Fisica Atomica, 
Molecular y Nuclear and Departamento de Ingenieria Electronica and Institute de 
Microelectronica de Barcelona (IMB-CNM), University of Valencia and CSIC, Valencia, 
Spain 

Department of Physics, University of British Columbia, Vancouver BC, Canada 
Department of Physics and Astronomy, University of Victoria, Victoria BC, Canada 
Department of Physics, University of Warwick, Coventry, United Kingdom 
Waseda University, Tokyo, Japan 

Department of Particle Physics, The Weizmann Institute of Science, Rehovot, Israel 
Department of Physics, University of Wisconsin, Madison WI, United States of 
America 

Fakultat fiir Physik und Astronomic, Julius-Maximilians-Universitat, Wurzburg, 
Germany 

Fachbereich C Physik, Bergische Universitat Wuppertal, Wuppertal, Germany 
Department of Physics, Yale University, New Haven CT, United States of America 
178 Yerevan Physics Institute, Yerevan, Armenia 

Centre de Calcul de I’lnstitut National de Physique Nucleaire et de Physique des 
Particules (IN2P3), Villeurbanne, France 

“ Also at Department of Physics, King’s College London, London, United Kingdom 
^ Also at Institute of Physics, Azerbaijan Academy of Sciences, Baku, Azerbaijan 
Also at Novosibirsk State University, Novosibirsk, Russia 
Also at TRIUMF, Vancouver BC, Canada 

® Also at Department of Physics, California State University, Fresno CA, United States of 
America 

^ Also at Department of Physics, University of Fribourg, Fribourg, Switzerland 
^ Also at Tomsk State University, Tomsk, Russia 

^ Also at CPPM, Aix-Marseille Universite and CNRS/IN2P3, Marseille, France 
* Also at Universita di Napoli Parthenope, Napoli, Italy 
^ Also at Institute of Particle Physics (IPP), Canada 
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^ Also at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United 
Kingdom 

* Also at Department of Physics, St. Petersburg State Polytechnical University, St. 
Petersburg, Russia 

Also at Louisiana Tech University, Ruston LA, United States of America 
” Also at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona, Spain 
° Also at Department of Physics, National Tsing Hua University, Taiwan 
^ Also at Department of Physics, The University of Texas at Austin, Austin TX, United 
States of America 

^ Also at Institute of Theoretical Physics, Ilia State University, Tbilisi, Georgia 
^ Also at CERN, Geneva, Switzerland 
Also at Georgian Technical University (GTU),Tbilisi, Georgia 

* Also at Ochadai Academic Production, Ochanomizu University, Tokyo, Japan 
“ Also at Manhattan College, New York NY, United States of America 

Also at Institute of Physics, Academia Sinica, Taipei, Taiwan 
Also at LAL, Universite Paris-Sud and CNRS/IN2P3, Orsay, France 
^ Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, 

Taipei, Taiwan 

y Also at Laboratoire de Physique Nucleaire et de Hautes Energies, UPMC and 
Universite Paris-Diderot and CNRS/IN2P3, Paris, France 

^ Also at School of Physical Sciences, National Institute of Science Education and 
Research, Bhubaneswar, India 

““ Also at Dipartimento di Fisica, Sapienza Universita di Roma, Roma, Italy 

Also at Moscow Institute of Physics and Technology State University, Dolgoprudny, 
Russia 

“ Also at Section de Physique, Universite de Geneve, Geneva, Switzerland 
Also at International School for Advanced Studies (SISSA), Trieste, Italy 
Also at Department of Physics and Astronomy, University of South Carolina, 

Columbia SC, United States of America 

Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China 
Also at Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, Russia 
Also at National Research Nuclear University MEPhI, Moscow, Russia 
Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, 
Budapest, Hungary 

Also at Department of Physics, Oxford University, Oxford, United Kingdom 
Also at Institut fiir Experimentalphysik, Universitat Hamburg, Hamburg, Germany 
Also at Department of Physics, The University of Michigan, Ann Arbor MI, United 
States of America 

Also at Discipline of Physics, University of KwaZulu-Natal, Durban, South Africa 
Also at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia 

* Deceased 
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